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Abstract 

The energy spectrum of high-energy neutrinos reported by the IceCube collaboration shows a 
dip between 400 TeV and 1 PeV. One intriguing explanation is that high-energy neutrinos scatter 
with the cosmic nentrino background through a ~ MeV mediator. Taking the density matrix 
approach, we develop a formalism to study the propagation of PeV nentrinos in the presence of 
the new neutrino interaction. If the interaction is flavored such as the gauged model 

we consider, the resonant collision may not suppress the PeV neutrino flux completely. The new 
force mediator may also contribute to the number of effectively massless degrees of freedom in the 
early universe and change the diffusion time of neutrinos from the supernova core. Astrophysical 
observations such as Big Bang Nucleosynthesis and supernova cooling provide an interesting test 
for the explanation. 
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I. INTRODUCTION 


The IceCube experiment has recently reported the observation of neutrinos in the energy 
range of TeV-PeV which provides the first evidence for extraterrestrial high-energy 
neutrinos. An interesting feature of the observed spectrum is a null detection of high- 
energy neutrinos in the energy range of 400-800 TeV. Although at present statistics have 
not been sufficient enough to confirm the existence of the dip in the spectrum, there have 
been investigations whether it can be explained by some new physics j5|-16|. One possibility 


is that the high-energy neutrinos may scatter with the cosmic neutrino background (Ci/B) 
and lose their energy, resulting in the dip-like feature in the spectrum 

This scenario has several interesting implications. To suppress the neutrino flux in the 
range of 400-800 TeV as indicated by the IceCube observation, the scattering cross sec¬ 
tion between the high-energy and Cz/B neutrinos must be signihcantly large, which can be 
achieved by the Breit-Wigner resonance. Since the resonance mass mR is close to the center- 
of-mass energy, it can be estimated as ttir ~ ~ 0(1) MeV, where neutrino energy 

El, ~ PeV and neutrino mass ruj, ~ 0.01 eV. Therefore, it predicts a new interaction in 
the neutrino sector with a force mediator much lighter than the weak scale. Because the 
resonant cross section is sensitive to Ei, and m,^, the IceCube neutrino ^ectrum may contain 
rich information about neutrino mass and redshift of the source 11. lil . 


From the perspective of particle physics model building, it is quite challenging to extend 
the lepton sector of the Standard Model (SM) with an additional interaction. For example, 
if the light mediator couples to the three generations of leptons universally, there are str ong 
constraints on the interaction strength from such as electron beam-dump experiments [171] 
and rare decays of mesons [I^. Therefore, it is reasonable to consider models in which the 
new interaction is not flavor-blind. In this case, it is important to treat the propagation of 
the high-energy neutrinos properly, in order to calculate the neutrino flux at IceCube. In this 
paper, we use the density matrix approach to study propagation of PeV neutrinos from the 
source to the IceCube detector in the presence of a new flavored neutrino self-interaction. To 
illustrate our main point, we consider an extension of the SM with a gauged — Tr Hi®- 
This model has several attractive features: it is gauge anomaly-free; it explains the nearly 
maximum mixing angle between the second and third generations; the model also evades 
severe constraints from electron beam-dump experiments. 

We also study cosmological and astrophysical implications of the model. Since the me¬ 
diator mass is close to the temperature of the Big Bang Nucleosynthesis (BBN) era, the 
presence of the mediator in the early universe can potentially contribute to the number 
of effectively massless degrees of freedom. In additional, the light mediator may also be 
produced in the core of supernovae. The frequent collision between neutrinos mediated by 
the new force may significantly reduce the neutrino mean free path, which slows down the 
supernova cooling process. We show that both BBN and supernova constraints are sensi¬ 
tive to the parameter region of the model explaining the dip of the IceCube PeV neutrino 
spectrum. Our result can be generalized to other models with an 0{1) MeV force carrier 
coupled to SM neutrinos. 

The rest of the paper is organized as follows. In the next section, we discuss the generic 
feature of gauged model and experimental constraints. In Sec. lIIIl we derive the 

Boltzma nn equation governing the evolution of the neutrino density matrix. Then, we dis¬ 
cuss cosmological and astrophysical implications in Sec.HVl Sec.|V]is devoted to summarizing 
our results. In APPENDIX, we present the derivation of the resonant scattering rate. 
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II. PARTICLE PHYSICS MODEL 


We assume SM leptons have a new interaction with the following Lagrangian 


= -g'z'^Y. [La^^ieLe + lR^^weR] , 


( 1 ) 


where g' is the coupling constant, Z'^ is the new gauge boson, Li denotes the lepton doublet, 
Ir denotes the lepton singlet, and the charge matrix is = diag(0,1, —1) in the interaction 
basis I = (e, /i, r). We furthermore assume that the gauge boson mass is mz' ~ MeV. 

There are several experimental constraints on this model. The existence of the light Z'^ 
opens up new decay channels for W and Z bosons such as three-body decays —?• h'^pZ'. 

These new processes chan ge W/Z-boson decay branching ratios by AT/T = in 

the limit of mz' <C my/jz [28|. To achieve ~ 1% precision of measured W/Z-boson decay 
branching ratios 0 , we estimate g' < 0.7. 

One of the most stringent constraints on g' is from the precise measurement of muon 
anomalous magnetic moment Uj, = {g — 2)/2. The leading contribution from Z'-exchange to 


can be evaluated as [19 


26,^ 


= 


V2 


Stt^ 


dx 


2m‘tx‘^{l — x) 


x‘^m? + (1 


xjrrir. 


( 2 ) 


In fact, several experiments have reported the measured value deviates from the SM predic 


tion at the level of Aa'^^^ = (42.6 ± 16.5) x 10 [mI, |32|. While the hadronic uncertainty is 


still in debate, we require Aa^ to be less than Aa^^'^, which gives rise to < 5 x 10“"^. 

The measurements of neutrino-electron interactions also put stringent constraints on 
the model. Although the Z' boson does not couple to electrons directly in our model, it 
contributes to neutrino-electron scattering through photon-Z' mixing radiatively induced by 
vacuum polarization with /r and r in the loop. Ref. jssj analyzed data from Borexino (solar 
neutrino) and GEMMA (reactor neutrino) jssj experiments to put constraint on the 
gauge coupling constant gs-L of a gauged B — L model, which is flavor-blind. To apply their 
result to our model, we hrst relax Borexino constraint on g-Q-h by a factor of (1/0.66) 
because the Z' boson couples only to and Vr, but not z/g, which accounts for about 34 % of 
the total solar neutrino flux [1^. We then impose the scaled upper bound on ^/eeg', where e 
is the electric charge of the electron, and e is the photon-Z' mixing parameter. We calculate 
the mixing parameter e as 


e = 


eg 


2tt^ 


dxx{l — x) In 


— x{l — x)q^ 
m'i — x{l — x)q^ 


( 3 ) 


where and m. 


are the masses of /i and r, respectively. We take a typical value of 
momentum transfer in neutrino-electron scattering, q^ = —IMeV^. Note the choice of q^ 
does not change e as long as \q^\ -C We hnd that constraint from Borexino experiment 
is more stringent than that from muon anomalous magnetic moment in the light Z' boson 
region of mz' ^ 10 MeV. Since the GEMMA experiment looks for reactor z/g’s before they 
oscillate (the distance from reactor is 13.9 m), it is not applicable to our model because Z' 
does not couple to z/g. 

The realization of the observed neutrino masses and mixing angles in the gauged — Lr 
model has been discussed in the literature 23-27|. In this paper, we assume that neu¬ 
trino masses are quasi-degenerate, which can be achieved with a proper choice of the 
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symmetry breaki ng pat tern j26[|. In this case, we can translate the cosmological limit 


my. < 0.25 eV 3^38] to an npper bonnd on the individnal nentrino mass m^. < 0.083 eV. 
On the other hand, the observed atmospherical nentrino mass is ~ 0.048 eV, 

which leads to a lower bonnd on m^. ~ 0.05 eV for a degenerate nentrino mass spec- 
trnm. Since the dip of the PeV nentrino spectrnm at IceCnbe is in the energy range of 
4 X 102 < Ey < 8 X 102 TeV 0, 8 and the resonance condition for nentrino scattering is 
m|, ~ 2E^miy^, we obtain a preferred range of Z' mass 5 MeV < mz' < 10 MeV. If scattering 
occnrs at a high redshift, we can shift this mass range by a factor of vTT^ accordingly. 


III. NEUTRINO PROPAGATION 


We hrst estimate the coherence length of PeV nentrinos as follows [39l-l41|. For any two of 
the nentrino mass eigenstates composing a flavor state, the velocity difference of their wave 
packets is \vi — Vj\ ~ |Am2.|/2ii^2^ where = rrif — m^. After they travel distance L, the 


wave packets are ~ L| Am 2 ^|/ 2 ii ^2 apart. If L|Am 2 ^|/ 2 ii ^2 jg larger than the nncertainty in 
their spatial location, the wave packets do not overlap and lose coherence^ Therefore, the 
coherence length can be estimated as 


^COh,2ji - 


|Am2^.| 




( 4 ) 


where is the spatial nncertainty of PeV nentrinos. Assnming that IceCnbe PeV nentrinos 
are prodnced by decays of high-energy pions, we expect the spatial nncertainty of PeV 
nentrinos is of the order of the distance that the pion travels before it decays, he., ax — 
/[AE^) [ 4 ^, where ttIt^/^AE^) is the Lorentz contraction factor, and the pion lifetime 
and mass in the rest frame are r.^- ~ 2.6 x 10“®s, and m^^ — 140 MeV, respectively. Taking 
E^=l PeV and |Am2^.| = 10“^ eV2, we can estimate the spatial nncertainty ax — 2.7 x 10“® 
cm and the coherence length Lcoh,p — 100 Gpc. Therefore, if PeV nentrinos are prodnced 
by decay of free pions, the coherent length can be larger than the particle horizon size of the 
nniverse ~ 14 Gpc. However, the environmental effects of the PeV nentrino sonrce, snch 
as collisions of the parent particle with particles in medinm and the presence of magnetic 


helds, may shorten the coherence length signihcantly |40l . I42l . I43l |. Fnrthermore, even if 


the coherent oscillation is maintained dnring the propagation, it may not be detected dne 
to nncertainties of sonrce distance and limitations of detector resolntion. In onr analysis, 
we hrst assnme that PeV nentrinos are coherent and derive the probability matrix in the 
presence of the new interaction, and then take time-averaging over the oscillatory terms to 
inclnde the possible decoherence effects. Therefore, the formalism we will develop below is 
valid even if PeV nentrinos have a coherent length shorter than the propagation distance. 

To describe propagation of the PeV nentrinos from the sonrce to the IceGnbe detector, 
we consider evolntion of the following density matrix 


oc {a\{k)ae{k )), 


( 5 ) 


where a^{k) is an annihilation operator of lepton nentrino with momentnm k. We normal¬ 
ize the density matrix snch that the nnmber density is given by ni{t) = J d^k/{27r)^^u{k, t). 


^ In this case, the density matrix defined in Eq. m is diagonal in the mass basis. 
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The evolution equation of the density matrix can be derived from nonequilibrium field the¬ 


ory 


44-46 


t) - t) = -i[J^{k), ^{k, t)] + ‘^[^], 

OT dk 


( 6 ) 


where H is the Hubble expansion rate, [•, •] denotes the commutator, Jif{k) is the Hamilto¬ 
nian, and represents the collision term. The Hamiltonian is 


J^{k) = ~ ijj{k) -I- {2u{k )), (7) 

where k = |fc|, u:{k) = \/k‘^ + m^, = tr(.^*^j,)/3, and — Since 

mass-squared differences are small, we have Am^ -C u}{k). 

We write the collision term as 
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( 8 ) 


and the total scattering rate is given by 


Vs{kA) = X{‘i)Ty{t)—'^5[k-m\,l{2my)\ , (9) 

rriy mz' 

where C(s) is Riemann zeta function. We present the derivation of Eq. ([9]) in APPENDIX. 

Taking redshift z and incoming momentum k^ as time and momentum coordinates instead 
of cosmic time t and physical momentum A;, respectively, we obtain the density matrix 
evolution of PeV neutrino. 


z) = z)] - z),M} 


where 


= <^(k,t), Aje = tffz, ^ ■ r. = r.(z)s 

2fco(l -h zYH[z) 


1 + z 


m 


Z' 


2myko 


and the optical depth Ts{z) is 


Ts{z) = 18C(3)Tj^o(l + z)^ 


1 P 


Z' 


H{z)mz^ mz' 


( 10 ) 

( 11 ) 

( 12 ) 


Integrating both sides of Eq. flTOj) from Zi to Zf, we obtain a formal solution for ^ 

^{ko, Zf) = ^{ko, Zf, Zi)^{ko, Zi)^(ko, Zf, Zi)^ (13) 

with non-unitary operator 

-i [ dz" 

Jz' 

where V is the propagation order operator dehned such that V {Q{z)Q'{z')} = Q{z)Q\z') 
for z < z' and Q'{z')Q{z) for z' < z. Substituting Eq. flTT]) into Eq. flTT]) . we obtain 

^{kQ,z',z) = ‘W{ko,z',Zs{ko))3"s{zs{ko),z',z)^{ko,Zs{ko),z ), (15) 


Jf{ko,z")-iTs{ko,z")^/2 


(14) 


^{ko, z', z)=V< exp 
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where 


^ (/co, z) = exp 




-i- 


2kn 


L{z',z) 


with L{z', z) = 


dz" 


(1 + z")m{z") 


, ( 16 ) 


^s{z'',z',z) 

Zs{ko) 


exp {-Ts{z'')M/2) (z' < A' < z) 

1 (otherwise) ’ 

= m|,/(2m^A;o) - 1 • 


(17) 


U is the unitary matrix that relates the mass basis z/j to a basis of interest vi = 'Yui 
the density matrix a.t z = Zf can be written as 





ArrA i ~ 


Amf ,■ ~ 

X exp 

-^-^L{zs,z,) 

AKq 

exp 



(18) 


For PeV neutrinos, the oscillation length Fosc.p = 4:Tiko/Amfj is 8x10 ®(fco/PeV)(eV^/m7) 

while the propagation length L{z', z) is on the order of l/iPo = 3 x lO^Mpc/h as long as 
A and z — z' are on the order of unity. This implies that we can take the period aver¬ 
age of the exponential terms, which gives rise to ^exp[—zZ\m^jT/(2fco)]^ — With this 
approximation, we obtain the density matrix as 


Fs,2^S 

(19) 


For PeV neutrino detection, it is useful to take the interaction basis. Here, the unitary 
matrix is called Maki-Nakagawa-Sakata (MNS) matrix ^ that is often parametrized by 


Vmns — 


■ 1 

0 

0 ■ 


0 

C23 

S23 


0 

— S23 

C23 



ci3 0 sise 
0 1 0 

-sise*'^ 0 Ci3 



Cl2 

Sl2 

o' 


-Sl2 

Cl2 

0 


0 

0 

0 


diag(l, ,(20) 


with Cij and Sij denote cos(6'jj) and sin(6*jj), respectively. In the interaction basis, the 
interaction matrix is diagonal. 




z ,z = 


( diag(l, 

0 


g-r42")/2^g-r.(^")/2) 


(z' < z" 


< ^) 
(otherwise) 


( 21 ) 


The elements of the probability matrix are 

= ^u{ko, Zf) for <^(^£ 200 , Zi) = 6 £^ei 6 £ 2 e'\{i^e{zf)\i'e'{zi))\'^-i (22) 

where we do not sum over £' and i. We can write explicitly as 

^ ^ Zf, Zi)U£^i-Jj£ij^Jj£ii-Jj£^j^^A/'s^i 2 i- 2 i.Zs-i Zf, Zi)U£^i^U£^^. (23) 

.£^S,2^S 


pc. 
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To evaluate ^ numerically, we take the following values for the MNS matrix parameters 
(normal hierarchy) [1^: sin^( 6 *i 2 ) = 0.308, sin^( 6 * 23 ) = 0.437, sin^( 6 *i 3 ) = 0.0234, and f^/vr = 
1.39. Combining Eqs. fl2l|) and ([23]), we obtain 


"0.30 0.13 0.12' 


0.07 -0.05 -0.03' 


'0.18 0.15 0.12' 

0.13 0.06 0.05 

_|_ ^-Ts{Zs)I2 

-0.05 0.03 0.02 


0.15 0.29 0.31 

0.12 0.05 0.04 


-0.03 0.02 0.01 


0.12 0.31 0.35 


(24) 


for Zs > 0 . 

With the probability matrix given in Eq. (I2T|) . we check several extreme cases. In the 
absence of scattering, i.e., Ts{zs) = 0, the flavor composition of PeV neutrinos at the IceCube 
detector is completely determined by the initial condition and oscillations. From Eq. (I2T|) . 
we can see that = 1; where i'{i) = e,/r, r. This is expected because without 

scattering, the total probability for Ending neutrinos in difierent flavors is conserved. 

In the limit of Ts{zs) ^ 1, one might think that and would be completely depleted 
in the neutrino flux reaching the IceCube detector because of collisions mediated by the Z'. 
However, this is not the case. Even though the last two terms of Eq. fl2T|) vanish, the first 
one does not depend on Ts{zs) at all. Therefore, the probability for Ending and Ur does 
not vanish. In general, for a model with flavor-dependent neutrino interactions, the CuB can 
not completely absorb each flavor of high-energy neutrinos. 

To check whether the model can produce the dip in the energy spectrum of 

high-energy neutrinos observed by the IceCube, we consider two possible sources for the 
high-energy neutrinos. If they originate from pp collisions, the initial flavor compositions 
in the neutrino flux are (0^^, 0^,^, 0^^) = ~ (1,2,0). Applying the probability 

matrix of Eq. fl24p with Ts{zs) ^ 1, we obtain the final flavor composition is (0.56, 0.25, 0.22), 
and the total survival rate is 1/3. While, for a py source, the initial flavor compositions 
are ( 1 , 1 , 0 ) and ( 0 , 1 , 0 ) for neutrinos and anti-neutrino, respectively. At the detector, they 
become (0.43,0.19,0.17) with a 2/5 survival rate and (0.13,0.06,0.05) with a 1/4 survival 
rate, respectively. Such a fractional suppression can explain the dip in the high-energy 
neutrino spectrum observed in the IceCube detector (see Fig. d). Within the current data 
set, the expected number of events is just 2. Therefore, even a 2/5 suppression factor can 
reduce the expected number of events below the Poisson limit, 1. On the other hand, our 
result can be tested with better statistics in the accumulated data set in the near future. 

To examine the favored parameter region for the — Lr model, we use Eq. ca and 
take mz' = 8 MeV, niu = 0. 05 eV, and ko = 600 TeV. The Hubble expansion rate is 
H{z) = + zY + Oa, where Hq = lOOh km/s/Mpc, VLm = 0.32, Oa = 0.68, and 

h = 0.69. We End the redshift at which the collision occurs is Zg — 0.07, and the optical 


depth is 


r, ~ 1 


9 


1.7 X 10 


-4 


(25) 


We see that the required value g' for the model to explain the dip is below the constraint 
from the muon g — 2 measurement, g' ^ 5 x 10“^. It is interesting to note that the model 
can explain the dip in the neutrino spectrum and the discrepancy in the muon anomalous 
magnetic moment. If we take g' = 5 x 10“^ as preferred by the optical depth is 

Tg — 8.5, which is more than enough to suppress the neutrino flux in the 400-800 TeV 
energy range at the IceCube detector [l^. However, as we will show later, this parameter 
region is strongly disfavored by the constraint from supernova cooling. 
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FIG. 1. Expected number of PeV neutrino events at the IceCube detector for the standard model 
(solid), the flavored model (dotted), and a model in which neutrinos have a new flavor- 

blind interaction (dashed), together with experimental data For the Lfj, — Lr model, we assume 
that PeV neutrinos originate from pp collisions. Because of the coherent effect in PeV neutrino 
propagation, the neutrino flux does not vanish near the resonance for the Eavored model. 


Our result can be easily generalized to the case in which the neutrino interaction is flavor- 
blind. After replacing ^ = diag(0,1,1) by diag(l, 1,1), we can write the probability matrix 
as 




0.54 0.24 0.21 
0.24 0.38 0.38 
0.21 0.38 0.41 


(26) 


In this case, all flavor compositions in the neutrino flux will be suppressed if > 1. 


IV. COSMOLOGICAL AND ASTROPHYSICAL IMPLICATIONS 


In this section, we study cosmological and astrophysical implications if neutrinos interact 
with a new ~ MeV force carrier. We take the L^ — Lt model as an example, and our analysis 
can be generalized to other models. 
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A. AAeff constraints 


In the early universe, Z' bosons can be produced in the SM thermal both by inverse 
decay and pair annihilation of leptons. The rate of inverse decay can be estimated as Tiny ~ 
X where mz' jT is the time dilation factor. The rate for pair annihilation is 

hann ~ ■ At early stages when the temperature is high, the pair annihilation process 

dominates the production of Z' in the thermal bath. But inverse decay becomes more 
important when T < 100 GeV(10“^/5'')(m2//10 MeV). When the temperature drops below 
m^/, the number density of Z' becomes suppressed by the Boltzmann factor. Since the Z' 
mass is close to the temperature when BBN starts, it may contribute to the effective number 
of neutrinos Wff. 

The mediator Z' may change Ws in two ways. If Z' is still relativistic during BBN, it 
contributes to Aefj directly. In this case, = 3 x 4/7 ~ 1.7, which is strongly disfavored 


48 


and CMB anisotropies j36l-l38| . Since 


by observations of light nuclei abundances 
mz’ ~ 10 MeV in the model we consider, the direct contribution to Agfr at T=0.1-l MeV is 
negligible because of the Boltzmann suppression factor. However, even in this case, Z' may 
still contribution to in an indirect way. When Z' becomes nonrelativistic, it transfers 
its entropy to and and increase their temperature relative to the temperature of v^. 
after neutrinos decouple from the SM thermal bath at Ti,^dec=l-5 MeV. To study this subtle 
effect, we take the following steps. We assume all neutrinos and anti-neutrinos have the 
same temperature Tj,dec=l-5 MeV when they decouple from the SM thermal bath. After 
decoupling, the and z/,-. 


and Z' form a thermal bath, which evolve independently from 


z/g and the photon. Then, we follow the phase space distribution function of z/^, z/,-, and 
Z' from Tj,dec=l-5 MeV to Ty =0.1 MeV, and derive a lower bound on mz' by demanding 
ANef[ < 0.7 at T =0.1 MeV. 

Since both inverse decay and pair annihilation processes respect CP, the relevant phase 
space distribution functions (per spin degrees of freedom) are given by 

1 1 

U = U, = fur = fur = fz' = . (27) 

where T' denotes the temperature of Vr and Z' after they decouple from the SM thermal 
bath, ^ and are the chemical potential per unit temperature for the neutrinos and Z' 
bosons, respectively. 

To evaluate ^ and we impose the following three conditions. 

• = 2^, because the inverse decay process is in the thermal equilibrium. 


The entropy per comoving volume is conserved, 

-f P P Sp^ + Sz')a^ = constant, 
where a is the scale factor. 


(28) 


The third condition depends on whether the pair annihilation process is in chemical 
equilibrium when the number density of Z' becomes negligible. We estimate the 
equilibrium condition requires g' > 10“®. In this case, Combining with the 

first condition, we have ^ = 0. If not, we instead use the conservation condition 
of the comoving number density 


(rz-i,^ + rip^ + + Up^ + 2nz')a^ = constant. 


(29) 
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We will discuss both cases. 


With the distribution functions given in Eq. fl27|) . we can write the entropy densities as 






(27r) 


sz' — 3 


Ank'^dk 

(27r)3 




m 


Z' 


V 


+ 


' Ak 

W' 




ST'^^k^ + ml, 


-r 


fz' 


(30) 

(31) 


In the case of = 0, we use the entropy conservation condition Eq. fl25]) to determine 

the temperature, T', at T=0.1 MeV for a given mz'- As we know, in the standard case, all 
neutrino species have the same temperature 0.1 x MeV when T=0.1 MeV. In our 

case, the Ue evolves as before, but both and Ur should have a higher temperature than 
0.1 X (4/11)3/3 because they inherit the energy density of the Z' boson. To evaluate the 
energy densities of the z/^, Ur and Z', we use 


Pl'ti Pvfj, Pvr Pvr 


A'lrk'^dk 

(27r)3 


kfu, Pz' — 3 


Ank'^dk 

(27r)3 


yV+m|i/z' ■ 


(32) 


We hnd that pz> is negligible mz' ^1 MeV at T=0.1 MeV. Using the standard dehnition. 


P'y + Pi/e + PPe + Pl^i, + Pv^ + Pvr + PPr 


P'1 


7 / 4 \ ^/^ 

! + -(-) (41V,,+ 3) 


(33) 


we calculate ANeg. Since the presence of the Z' boson does not change the thermal history 
of 7 and the following relation is still valid in our model 


P'y V Pue V Pl^e 


P'y 




(34) 


We demand dehned in Eq. fl33|) to be less than 0.7, and derive an upper bound m^/ > 
5.3 MeV(Ty dec/l-5MeV) shown in Fig. [2l We note that AN^^ drops signihcantly for larger 
rriz'- For example, AN^s is 0.1 for rriz' ^ 10 MeV(Tj, dec/l-5MeV). This is because the energy 
density carried by the Z' boson is suppressed by the Boltzmann factor ~ exp(—m^//T'). 

If the process of pair annihilation and creation becomes less than the Hubble expansion 
rate when the Z' boson becomes nonrelativistic, T', ^ and can be determined by solving 
Eqs. (|2^ and (12^ simultaneously with the initial condition ^ = ^' = 0 at T^^^dec = 1.5 MeV, 
where the number densities are given by 


— 


Aiik'^dk 

(27r)3 


/i/, nz' — 3 


Ank^dk 

(27r)3 


fz' 


(36) 


Following a similar procedure, we obtain mz' > 5.3 MeV(Tj,_dec/l-5MeV) for AN^s < 0.7, 
which is similar to the upper bound for the case of Therefore, our upper bound on 

mz> shown in Fig. |2] changes only a few percent even in the small coupling region. 
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B. Supernova Cooling and Neutrino Bursts 


The presence of a new MeV force carrier between neutrinos also has interesting implica¬ 
tions for the physics of supernova neutrinos. We hrst briefly summarize the basic picture 
in the standard case, see Refs. 50, IM| for review and references therein. The core-collapse 
supernova forms a proto-neutron star in its core. Its size and temperature are i? ~ 10 km 
and T ~ 30 MeV, respectively. In the core, nuclear reactions and electron pair-annihilations 
produce large numbers of neutrinos. These neutrinos reach thermal equilibrium with nuclear 
matter and cannot escape from the core due to its high density. As density and tempera¬ 
ture decrease with distance from the core, the mean free path of neutrinos becomes longer. 
Above some radius (called the neutrino sphere), they start streaming freely. The radius 
and temperature of the and sphere are roughly i? ~ 15 km and T ~ 8 MeV, respec¬ 
tively. The z/g sphere has a lager size and lower temperature since they can interact with 
circumstellar media more strongly through the charged current. The diffusion time can be 
estimated as Tdifr = A/c(i?/A)^, where A is the neutrino mean free path, and c is the speed 
of light. If neutrinos have only the SM weak interaction, we can estimate ^ 10 s, which 


is consistent with observed duration of the neutrino burst from SN1987A [52,153 


If neutrinos have new interactions, the standard picture of supernova neutrinos changes. 
For the L^—L^ model we consider, the Z' mediator can be produced inside the core if its mass 
is comparable or less than the core temperature. The the travel distance before it decays is 
only cTz' ~ 10“®km(5f'/10“^)“^(T/10MeV)(10MeV/m^/)^, which is much smaller than the 
core radius. Therefore, the produced Z' boson will reach thermal equilibrium with neutrinos 
and other particles in the core. These reactions may prevent Vy and Vj. from free-streaming, 
and change their diffusion time. We evaluate the diffusion time in the following way. The 
number density of and Vr increases as Z' decays, hy ~ l/2(r^/)nz/, where (T^/) is the 
total decay width of Z' averaged with phase space distributions. Meanwhile, riy decreases 
through the inverse decay. Therefore, we have hy = ll2(rz')nzi — {av)-^^^nyny. Since 
fij, = 0 in equilibrium, the detailed balance tells us {av) = n^,/(2 (T^/) where 

are the number densities when the particles are in thermal equilibrium as given in Eq. 0351) . 
Therefore, we can evaluate the mean-free path for z/„ and Vj. as 


A = 


{ll2{Vz')n^,ln7y 


(36) 


To estimate the mean free path in Eq. (]36|) . we factorize the thermally-averaged total width 
as {Tz') = g''^/{l2TT)mz'foimz'/T), where 


fo = Ibvr / Yi{k)fz'{k) / n(^ (1 - fy{p)) / n(p ) (1 - fy{p')) (27r) - p - p) (37) 


with n(fc) = d^k/2u;{27iy as the phase space measure. We evaluate fo numerically as shown 
in Fig.O Roughly speaking, fo can be regarded as the time dilation factor fo ~ mz'/T. 

Using Eq. 036p . we estimate the diffusion time for z/^ and z/^ in the presence of Z' as shown 
in Fig. |2] (solid black), where we have taken the core size as 10 km and temperature 8 MeV. 
For ruz' ~5-10 MeV and g' > 10“®, z/^ and Uj. may not contribute to the neutrino cooling 
of the core, and the neutrino burst would last 3 times longer than expected in the standard 
case. This appears to be incompatible with the observed ~ 10 s duration of the neutrino 
burst of SN1987A, although we can not draw a concrete conclusion because uncertainties in 
supernova modeling (e.g., nuclear equation of state), limited statistics of observed events. 
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FIG. 2. Summary of the parameter space of the gauged model. Experimental bounds on 

the model are from the muon g — 2 measurement (red band) and the electron-neutrino interaction 
measurement (brown dotted) at Borexino experiment (Sec.\I^. In the case of the quasi-degenerate 
neutrino mass spectrum, lower and upper bounds (blue vertical) on neutrino masses indicate the 
range of gauge boson mass mz' = 5-10 MeV required to produce the IceCube dip via resonant inter¬ 
action with CvB neutrinos fSec.tW). The optical depth (green) of the IceCube neutrinos fs should 
be larger than unity to reproduce the observed PeV neutrino spectrum fSe dllB) . Depending on 
parameters, resonant scattering can occur at different redshifts. Cosmological constraints requires 
mz' ^ 5 MeV such that < 0.7 (gray vertical). For mz' in the range of 5-10 MeV, the energy 

density carried by the Z' boson may still give sizable contributions to AN^ff, i.e., AN^g ^0.1-0.7 
fSec. lfV'31) . The resonant interaction may change the diffusion time of and (black). If 
is larger than ~ 10 s, it may delay supernova cooling (Sec. IIVBI) . 


and also uncertainties in deriving our limit. On the other hand, cooling processes through 
other invisible particles can compensate the su ppresse d neutrino cooling in this model. For 
example, QCD axions with mass of 
parameter of 10-10 


meV |54l-l56j and hidden photons with a mixing 


l57l. I58| are well-motivated candidates in charge of invisible cooling. 


Running simulations with this model are warranted for improving the limit and comparing 
it with observations in detailH 


^ The simulations are performed for QCD axions. 
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FIG. 3. Prefactor fo of the thermally-averaged decay width, {Tz') = /{12TT)mz' fniruz'/T), as 

a function of the ratio of Z' mass to temperature. By using this result, we evaluate the mean free 
path of Ufj, and Vr in the supernova core (see Eqs. if36)) and [37(\ ) . 


V. CONCLUSIONS 


The IceCube experiment has observed high-energy cosmic neutrinos for the hrst time. 
The observed neutrino spectrum exhibits a dip around the sub-PeV energy scale, which may 
indicate new physics beyond the SM of particle physics. One possible explanation is that 
PeV neutrinos may scatter with the Cz/B through a MeV resonance and lose their energy 
before reaching the IceCube detector. In this paper, we have developed a formalism to trace 
the propagation of PeV neutrinos in the presence of the new interaction. For the flavored 
interaction, we have shown that resonant scattering may not suppress the PeV neutrino flux 
completely, which could be tested in the near future. 

We have also discussed astrophysical and cosmological constraints on this type of models. 
The MeV mediator could be produced in the core of supernova, and frequent neutrino 
collisions induced by the mediator in the core could trap neutrinos inside the core. In the 
early universe, the mediator could thermalize with the SM thermal bath and contribute to 
the number of effectively massless degrees of freedom. We have shown both the BBN and 
supernova observations are sensitive to the favored parameter region explaining the dip in 
the IceCube high-energy neutrino spectrum. 
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APPENDIX 


The collision term describes scattering between the PeV neutrinos with Ci/B neu¬ 
trinos, I'lik) -|- VBV'i.'P) ^ which is given by 

[ dU{k’) [ dU{p) [ dU{^){27r)^6\k + p-k'-p')J2 


2 u{k) 
1 


spins 


^ 2 E *)• <)).<*}/»(?. t) (1 - /(.'(p". *)) 

^ b,b' 

-{1 - s^(k, t),^pk^{k', «).<*} (1 - Mp, t)) t) 


(38) 


where n(/c) = d^k/{2u{k){27r)^) is the phase space measure, {•, •} denotes the anticommuta¬ 
tor, is scattering amplitude. In deriving Eq. fl38|l . we have taken a basis in which 

the density matrix of Cz/B neutrinos is diagonal, {k,t) = 5^^' fb{k,t). 

The collision term '^[^] can be further simplified. Since the phase density of the PeV 
neutrinos is much less than the quantum limit, the Pauli blocking effect is negligible, he., 
1 — ^{k', t) ~ 1. We also neglect the inverse scattering process. In addition, we assume the 
distribution of background neutrinos are flavor-blind, he., fb{k,t) = f{k,t), where f{k,t) = 
l/[exp{k/Ti,(t)) + 1] with Ty(t) = -1- z) ~ 1.7 x -|- z) eV and z parametrizing 

the redshift. With these considerations, the collision term can be written as 

X E *), ^P'.^p}f(.p, t). (39) 

spins b,b' 


j dn{p){2TTf5\k + p-k'-p') 


In general, all s, t and w-channel exchanges of Z' contribute to neutrino scattering. Here, we 
focus on the parameter region where the s-channel resonance has a dominant contribution. 
In this case, the invariant amplitude matrix is 

= M {iy{k) -h P{p) u{k') + u{p')) , (40) 

with = ^bi^b'i' (see Eq. ([T])). Hence, we have 


^ = 2|A4 {u{k) + P(p) ^ u{k') + u{p')) , (41) 

b,b' 


where Note that M contains all information about flavor structure of new 

neutrino interactions. For the L^ — L^ model we consider, M = diag(0,1,1) in the interaction 
basis. Noting the definition of scattering cross section 

2 uj{k)2u{p)a{u(k) -|- v{p) v + h)nrei = j dll{k') J dIl(p)(2TrY5'^{k + p — k' — p') 

X ^ |A4 {i>{k) + Pip) iy{k') + P{p')) P , (42) 

spins 

we write the scattering rate Ts{k,t) as 

Ts{k,t) = 2 J t)a{iy{k) + P{p) ^ u + u)vpei ■ (43) 
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The cross section for resonant scattering is given by the Breit-Wigner formula 


ctr = An 


2J + 1 


(2si + 1)(2 s2 + 1) 


Bn^Br 


out 




t'cm 




<2 

cm 


m 


R. 


|2 + ELr^^cm) 


( 44 ) 


where si and S 2 are the spins of initial particles, J is the spin of the resonance, rriR is its 
mass, T/^ is its decay width, and and Bront are decay branching ratios to initial and 
hnal state particles, respectively. In the limit of ^ Eq. fl44)) can be written as 


an ~ 16n^ 


2J+ 1 


(2si + 1)(2 s2 + 1) 


Bn 


iBrout—5(^cm 

rriR 


(45) 


For the model we consider, Si = S 2 = 0, Brjn = Brout = 1/2, and we have 


a = 12w^^S(EL-’nl'), ( 46 ) 

mz' 

where Tz' = g'‘^mz'/{12n) is the Z' decay width in the rest frame. Therefore, the total 
scattering rate is 


r,(fc,t) = 9C(3)T3(t)^I^5[fc-m|,/(2m,)] , 
where ((s) is Riemann zeta function. 


(47) 
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